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摘 要 : 0. 618 法 是 一 维 线 搜 索 中 针对 一 维 单 峰 函 数 ， 应 用 最 为 广泛 的 一 种 方法 。 具 有 
良好 的 收敛 性 ， 但 其 收 敏 性 太 慢 ， 因 此 ， 本 文 基于 函数 在 搜索 区 间 端 点 和 区 间 内 任 一 点 函 
数值 的 基础 上 ， 给 出 了 一 种 普 适 性 的 线 搜索 加 速 策略 ， 每 步 选 代 都 可 以 在 较 大 程度 上 缩小 
子 数 值 的 不 确定 性 区 间 。 数 值 试 验 结果 表明 ， 其 收敛 速度 较 0. 618 法 有 所 提高 ， 尤 其 是 当初 
始 区 间 两 端 函数 值 相差 较 大 或 很 大 的 情况 下 ， 本 文 改进 算法 可 以 很 大 程度 上 减 小 区 间 范 围 。 
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An improved interval interpolation method based on 0.618 


method 


ZONG Yu-han*, ZHANG Yan-bo 
(School of Science,Beijing University of Posts and Telecommunications,Beijing 
100876,China) 


Abstract: The 0.618 method is the most widely used method for one-dimensional unimodal 
function in one-dimensional line search.It has good convergence, but its convergence 1s too slow. 
Therefore, based on the search of function values at the end points of the interval and any point in 
the interval, this paper presents a universal line search acceleration strategy. Each iteration can 
reduce the uncertainty interval of function values to a large extent.Numerical results show that the 
convergence speed of the proposed algorithm is better than that of the 0.618 method, especially 
when the function values at both ends of the initial interval differ greatly or greatly, the improved 
algorithm in this paper can reduce the interval range to a large extent. 
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strategy 


0 引言 然后 采用 某 种 插值 或 分 割 技术 缩小 这 个 区 间 ， 进 行 
精确 求解 。 此 类 方法 包括 黄金 分 割 法 、 分 数 法 及 二 次 
在 最 优化 的 迭代 算法 中 ， 需 要 用 到 一 维 搜索 或 。 ”插值 等 。 在 一 维 搜索 中 ， 一 般 要 求 函数 是 单 峰 函数 ， 
简称 线 搜索 ， 这 是 许多 非 线性 规划 算法 不 可 或 缺 的 因此 本 文 只 考虑 单 峰 函 数 情况 。 


一 部 分 。 它 可 归结 为 单 变量 函数 的 极 小 化 问题 ， 现实 中 大 多 数 函 数 的 导数 是 不 可 用 或 极 难 求 得 
其 数学 模型 为 : min f(x) 。 其 分 为 精确 线 搜索 和 非 ”的 ， 这 时 候 直接 法 的 应 用 便 广 泛 起 来 。 从 上 述 直接 
精确 线 搜索 于 ， 而 精确 线 搜索 又 分 为 直接 法 和 解 ”法 的 思想 可 以 看 出 ， 决 定 直接 法 性 能 的 关键 有 两 点 : 
析 法 。 解 析 法 是 使 用 函数 导数 的 搜索 ， 如 插值 法 、 牛 ” 一 是 初始 区 间 的 选取 ， 二 是 包含 极 小 值 点 子 区 间 的 
顿 法 及 抛物 线 法 等 ， 但 当 函 数 的 导数 不 可 用 时 ， 就 ”缩小 中。 

需要 采用 不 需要 导数 的 方法 ， 即 直接 法 求解 ， 其 基 0. 618 法 是 一 维 直 接 搜索 中 应 用 最 为 广泛 的 一 


本 思想 是 : 首先 确定 包括 问题 最 优 解 的 搜索 空间 ， 种 方法 ， 其 优点 是 对 函数 的 连续 性 不 做 要 求 ， 每 次 


友 代 只 计算 一 个 新 的 试探 点 ， 另 一 
具有 良好 的 收敛 性 ， 收 敛 速度 是 线性 的 。 但 缺点 是 
收敛 太 慢 。 本 文 在 0. 618 法 的 基础 之 上 ， 提 出 了 一 
种 改进 的 优化 方法 。 其 并 将 其 与 0. 618 法 比较 ， 通 
过 数值 计算 证 明了 此 种 方法 收敛 速度 快 ， 友 代步 数 


少 。 
1 0. 618 法 


0. 618 法 也 称 黄金 分 割 法 ， 其 基本 思想 四 是 
通过 对 每 步 试 探 点 函数 值 的 比较 ， 使 包含 极 小 值 点 
的 搜索 区 间 (不 确定 区 间 ) 不 断 缩 小 ， 当 达到 某 种 
准则 时 ， 可 以 找到 极 小 值 点 的 近似 。 

0. 618 法 取 试 探 点 的 规则 为 : 


1 =w +0.382(b, —a,) 


Li; = a +0.618(Db, 一 al) 
其 计算 步骤 如 下 : 


步骤 0: 给 定 初始 区 间 [46,b,] 及 精度 & <(0,]) 
计算 试探 点 和 局 ， 计 算 公式 为 
4 =a+0.382(D —a,) 
LW =a0o+0.618(b, 一 ao) 
计算 函数 值 f(4) 和 f(4),， 令 k=1。 
步骤 1: 若 满 足 |b. -ai|<s， 则 停止 计算 ,得 


到 极 小 值 点 的 近似 x = 


以 十 Qi 


否则 ， 当 


f(44)> Ho) 时 ， 转 步骤 2; 当 f(4)<f(p) 时 
转 步 又 3。 


步骤 2: 置 4i, = ， 


bin =b, Nin= Hi, 


Hi = Gil 二 0.618(D — AQ), 


计算 函数 值 FL )， 


转 步 又 4。 


步骤 3: 置 ar = 和，p =H Hen = hi, 


Ma = Aart+0.382(D, 一 ae) ， 


个 是 重复 利用 的 ; 


计算 函数 值 7( 4 ,) ， 转 步 又 4 
步骤 4， 置 上 := 大 +1， 返 回 步 骤 1 


2 基于 区 间 插 值 的 优化 算法 


文献 [6] [7] 提 出 了 基于 区 间 缩 小 的 不 同 改进 算 
法 ， 它 们 分 别 利 用 区 间 一 端 一 阶 导 数 和 二 阶 导 数 ， 
进一步 缩小 区 间 ， 但 对 于 导数 不 可 求 或 极 难 求 得 时 ， 
上 述 方法 就 不 再 适用 。 因 此 在 本 文中 ， 我 们 希望 通 
过 利用 当前 搜索 区 间 两 端点 的 函数 值 及 此 区 间 上 任 
意 一 点 的 函数 值 ， 构 造 出 一 种 一 维 线 搜索 的 加 速算 
法 ， 来 改进 0. 618 法 。 


设 f(x) 为 导数 不 可 求 或 极 难 求 得 的 一 维 连续 


单 峰 函数 ，[ai ,bi | 为 包含 极 小 值 点 x* 的 当前 搜索 


区 间 ， 首 先 我 们 选取 f(a ) 与 f(b, ) 中 函数 值 大 的 


点 ， 比 如 f(a )< f(b,)， 则 一 定 存 在 一 点 6 满足 


JE)= f(a), € ela,b,] 


上 述 属 性 表明 区 间 必 ,Bb | 不 包含 x* 。 所 以 我 们 可 以 


在 进行 下 一 步 前 如 先 去 掉 区 间 [E,b.]， 之 后 再 对 新 
， 会 使 算法 具有 更 好 的 效能 。 若 想 利 
用 上 述 优点 ， 关 键 在 于 如 何 求 使 得 (6)= f(a )。 
若 所 给 问题 函数 简单 ， 可 以 得 到 满足 条 件 的 精确 解 ， 


但 适用 范围 很 局 限 ， 此 时 就 需要 另 一 种 求解 5 的 通 
用 方法 。 
因为 插值 法 是 只 利用 函数 值 信息 ， 不 涉及 导数 


信息 ， 因 此 我 们 提出 用 插值 法 来 求 得 可 代替 & 的 一 
个 近似 值 5 。 具 体 思想 为 : 利用 函数 在 初始 区 间 
[40,5 |] 端点 处 的 函数 值 (ao)，f(6b6) 及 初始 区 间 


内 任 一 点 ce (a0,b,) (一 般 我 们 选 此 点 为 区 间 中 


点 ) 的 函数 值 f(c)， 来 构造 目标 函数 的 二 次 交 近 


函数 G(x)， 即 利用 三 点 二 


次 插值 方法 构造 遇 近 函 


数 四 ， 之 后 利用 此 函数 求解 满足 Wc)= 9(ai ) 的 6 。 
体 过程 加 为 ， 通过 曲线 = f(x) 的 三 点 


人, 月，(e, 万)，(e, 广 ) 作 二 次 函数 了 = p(x)， 


a)+a (f(a 
fb. Xa -ci)+ fla, Xe -bi)+ 


对 于 G(6)=G(B,)， 类 同 。 
下 面 给 出 此 改进 方法 的 人 


法 步骤 : 


步骤 0: 给 定 初始 搜索 区 间 [ao,b, | 和 区 间 内 一 
ao+bo 二 本 
记名 半 及 精度 se (0,1) 。 计 算 函 数值 (ao)， 


Jo)， oe 置 大 =0。 


步骤 1; 若 满 足 |b. -ai|<e， 则 停止 计算 ， 得 


已 Ta 


到 极 小 值 点 的 近似 x = 本 。 否则 转 步 骤 2 


步骤 2: 若 Fa)< f(b;)， 按 公式 


Ve 


/ee eee ay+f(c)a -bh) 


fb Jar -ci)+ flar ke 


(f(a Nbi-e 人 
-b,)+f(c, 和 


计算 G ， 去 掉 区 间 [s,b.]， 置 b=6 得 新 的 区 间 


[a,,5b, ]， 转 步骤 4。 否 则 转 步 又 3 


步骤 3: 当 f(a; )> f(b,) 时 ， 按 公式 


攻 (eoles -b+ (0, je 7 "+ /lo )- a)+h 


(fabi-e) + foe -a)+ fecha -b)) 
TO er ee 所-a) 


hj+ fl 


计算 9G ， 去 掉 区 间 [a4,6]， 置 a =5 得 新 的 区 间 


[c ,六 ]， 转 步骤 4 


步骤 4: 计算 试探 点 轴 和 yi， 公式 为 
1 = a, +0.382(b, —a,) 
Mi =a, +0.618(b, —a,) 

计算 函数 值 f() 和 f(44)。 当 了 (4)> Au) 时 ， 


转 步骤 5; 当 f(4)< (4) 时 ， 转 步骤 6 


步骤 5: 置 4, =，b1=b, c= 
转 步 又 7 

步骤 6: 置 4j, =a1,， Di=W, Cin =h， 
转 步 又 7 


步骤 7: 置 上 :=k+1， 返 回 步 又 1。 


3 数值 实验 

为 验证 改进 算法 的 效果 ， 我 们 选取 一 个 简单 的 
一 维 单 峰 函 数 和 两 个 较 复杂 的 单 峰 函数 外 来 比较 改 
进 的 算法 与 0.618 法 ， 将 它们 对 例 1、2、3 中 目标 函 
数 在 指定 区 间 进 行 搜索 ， 通 过 对 比 两 算法 的 迭代 步 


es 其 中 精度 
我 们 分 别 选 取 102 ，104 ， 结 果 如 表 2、 表 3、 表 4 所 
示 。 


例 1: f(x)=x ， 最 小 值 点 为 0， 初 始 区 间 选 为 [- 


1, 100]. 


例 2: f(x)=e' -3x ， 最 小 值 点 为 -0. 144275， 初 
台 区 间 选 为 [-1, 4]. 
例 3 : f(x)=cosx+(x-2), 


2. 35424， 初 始 区 间 选 为 [1 4]. 


最 小 值 点 为 


表 1 函数 1 在 两 种 算法 下 的 结果 比较 
Table 1 Comparison of the results of 
function 1 under the two algorithms 
改进 算法 
[-1 100] 


0.618 法 
[-1 100] 


初始 区 间 


迭代 步 数 6 
输出 结果 1.15431x103 
初始 区 间 [-1, 100] 
迭代 步 数 10 
输出 结果 -9.13146x10” 


20 
2.11656x103 
[-1, 100] 
29 
7.23312x105 


表 2 函数 2 在 两 种 算法 下 的 结果 比较 


Table 2 Comparison of the 


results of 


function 2 under the two algorithms 


0. 618 法 
[-1, 4 
10 
-0. 144297 
[-1, 4 
23 


改进 算法 
初始 区 间 [-1, 4] = 
迭代 步 数 5 
输出 结果 -0. 147635 
初始 区 间 [-1, 4] 
和 欠 代 步 数 13 
输出 结果 -0. 148459 


-0. 144314 


表 3 函数 3 在 两 种 算法 下 的 结果 比较 


Table 3 Comparison of the 


results of 


function 3 under the two algorithms 


改进 算法 
初始 区 间 [1, 4] 
友 代 步 数 6 
输出 结果 2. 35115 
初始 区 间 [1, 4] 
迭代 步 数 12 
输出 结果 2. 3468 


0.618 法 
[1, 4] 
11 
2. 3542 
[1, 4] 
22 
2. 3542 


通过 观察 三 个 函数 的 测试 结果 ， 改 进 算法 的 迭 
代步 数 丰 仅 比 0.618 法 的 娄 小 和 多， 而 且 最 优 值 点 
精度 也 有 明显 改进 ， 同 时 构造 出 的 点 列 能 够 在 有 限 


步 内 得 到 达到 所 要 求 精度 的 最 优化 问题 的 最 优 解 。 
因此 改进 算法 不 仅 有 稳定 的 收敛 性 


， 而 且 在 各 种 情 


况 下 也 能 保持 较 快 的 收敛 速度 。 


同时 经 过 分 析 可 知 ，0. 618 法 每 步 迭 代 只 计算 


一 个 函数 值 ， 而 改进 的 和信 
但 因为 迭代 步 数 减少 不 
有 一 定 提 


1 | 
alm i 示 


法 每 步 计 算 两 个 函数 值 ， 
因此 改进 算法 也 是 
高 的 。 我 们 通过 测试 函数 说 明 这 一 结论 ， 


针对 例 1 在 不 同 长 度 的 区 间 上 的 区 间 缩 短 率 ， 对 


0.618 法 


和 改进 算法 进行 比较 ， 结 果 如 表 4 所 示 。 


表 4 两 算法 对 于 不 同 区 间 的 计算 
Table 4 Two algorithms for different 


interval calculation 


精度 : 0.618 法 改进 算法 
10” 友 代 输出 迁 代 输出 
初始 步 数 结果 步 数 结果 
区 间 
[-1, 500] 33 -3.96339x105 10 -9.13188x10 
[-1, 100] 29 -7.23312x10% 10 -9.13146x10 
[-1, 10] 25 -2.76438x10% 10 -9.08866x10 
1; 如 21 -4.40305x105 ”10 -8.88847x10 
[-1, 0. 1] 20 -3.79423x105 8 -6.25798x10 
[-1, 0. 01] 20 -4.57181x105 5 -1.86693x10- 


观察 表 4， 当 初始 区 间 是 [-1 500] 时 ， 友 代步 
数 是 0. 618 法 的 三 分 之 一 倍 ， 同 时 区 间 缩 短 率 减 小 
了 28%; 当初 始 区 间 是 [-1, 0.01] 时 ， 迭 代步 数 是 
0.618 法 的 四 分 之 一 倍 ， 同 时 区 间 缩 短 率 减 小 了 
36%; 而 当初 始 区 间 是 [-1, I 时 ， 和 迭代 步 数 虽然 是 
0. 618 法 的 五 分 之 二 倍 ， 但 是 区 间 缩 短 率 只 减 小 了 
6%。 因 此 我 们 可 以 发 现 当 初始 区 间 两 端 函数 值 相差 
较 大 或 很 大 的 情况 下 ， 比 如 初始 区 间 两 端点 相差 
100 倍 左右 时 ， 改 进 的 算法 可 以 较 快 的 减 小 区 间 范 
围 ， 以 较 快 的 速度 收敛 到 最 小 值 点 。 

当 区 间 长 度 慢 慢 变 小 时 ，0. 618 法 的 迭代 步 数 
变化 很 小 ， 或 者 几乎 不 变 ， 但 是 相 比 较 而 言 ， 改 进 
的 算法 在 区 间 长 度 减 小 时 ， 不 会 出 现 停滞 的 现象 ， 
且 达 代 步 数 还 会 偶尔 成 倍 减 小 ， 因 此 改进 的 算法 可 
以 显现 出 很 好 的 性 能 。 


4 结语 


线 搜索 算法 是 单 因素 实验 设计 中 最 常用 的 一 种 
方法 ， 其 收敛 速度 直接 影响 算法 的 效率 。 如 何 利用 
函数 在 搜索 区 间 的 信息 ， 提 高 算法 的 收 合 速 度 是 相 
关 研 究 的 重点 。 本 文通 过 利用 目标 函数 在 搜索 区 间 
端点 处 的 函数 值 信息 ， 不 断 缩小 每 步 迭 代 区 间 ， 获 
得 了 具有 普 适 性 的 加 速 策 略 ， 进 行 的 数值 实验 表明 
本 研究 所 给 出 的 算法 能 够 加 快 现 有 线 搜索 算法 的 收 
敛 速度 。 
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